Critical power of collapsing vortices 
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We calculate the critical power for collapse of linearly-polarized phase vortices, and show that 
this expression is more accurate than previous results. Unlike the non-vortex case, deviations from 
radial symmetry do not increase the critical power for collapse, but rather lead to disintegration 
into collapsing non-vortex filaments. The cases of circular, radial and azimuthal polarizations are 
also considered. 
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. The nonlinear optical process of self- focusing sets an upper limit on the amount of laser power that can be propagated 
through a medium with an intensity dependent refractive index (i.e., n = no + TI2I, where no is the linear refractive 
, index, n2 is the nonlinear refractive index, and I is the intensity). For powers above this threshold the beam will 
undergo collapse, with the peak intensity becoming sufficiently high that damage to the material can occur. Ultimately, 
collapse will be arrested by some physical mechanism, such as plasma formation, normal dispersion or damping. 
Let us briefly review the situation in the non- vortex case. The value of the critical power is given by [l[ 
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where p cr is the non-dimensional critical power for collapse in the dimensionless NLS 

iipz(z,x,y) + Aip + \vj\ 2 ip = 0, ip(0,x,y) = i>o(x,y). (1) 

Oh| In the NLS model, there is no mechanism for arrest of collapse, hence collapse is defined as the maximal amplitude 
becoming infinite. Weinstein Q proved that the lower bound for the critical power is equal to p cr — j \R\ 2 rdr w 1.86, 
' i.e., the power of the Townes profile, which is the ground state solution of 
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R" + -R' -R+ R? = 0, R'(0) = 0, R(oo) = 0. 
>V r 



While the Townesian input beams ip = AiJ(Ar), where A > 0, can collapse with exactly the input power p cr , all other 
input profiles require power strictly above p cr for collapse [3, 4]. In practice, however, the critical power of peak-type 

(i.e., non ring-type) radially-symmetric input beams is only a few percents above p cr 0, ||. For example, the critical 

4 

power of Gaussian and super-Gaussian (-00 = ce~ r ) input beams is « 2% and as 8% above p cr , respectively. 

We now consider the critical power of vortex input beams. In Q, Kruglov et al. derived an expression for the 
j critical power of vortex beams, and showed that it increases with the winding number (or topological charge) m. 
^\ . In this study, we show that this expression is inaccurate, and derive the correct expression for the critical power. 
£NJ ' Unlike the vortex-free case, deviations from radial-symmetry do not increase the critical power, but rather lead to 
, disintegration into collapsing non-vortex filaments. 

■ We first consider radially-symmetric vortex input beams of the form tpo = Ao(r)e l . In this case, the solution 
remains a vortex with winding number m, i.e., it is of the form yb{z,r,6) — A{z,r)e lme Following a similar 
derivation to , it can be rigorously shown that the lower bound for the critical power of radially-symmetric vortex 
input beams tpo = A (r)e ml6 is 
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Pcr(m) = \R m \ rdr 



where R m is the ground state solution of 



1 / m 2 \ 

R' r ' n {r) + -R' m - (l + — j R rn + Rf n = 0, R' m (0) = 0, ^(00) = 0. 

The values of p cr (m) for m = 1, • • • , 6 are listed in Table |U Using the approximation 



R m (r) ^VSsecht - (2) 
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TABLE I: 



m 


1 


2 


3 


4 


5 


6 


Pcr(m)/p cr 


4.12 


7.65 


11.3 


15.0 


18.7 


22.4 




FIG. 1: (color online) Critical power p cr (m) = R^ n rdr (x), the approximation 4y/3m (dashed line), numerical estimate of 
the critical power as a function of m (o, data taken from [(|), the critical power for collapse of Laguerre-Gaussians (dash-dotted 
line) and the analytic estimate I^r (0i solid line). 



we can derive the analytic approximation p cr (to) « 4v3m. Figurc[I]shows that p cr {m) is well approximated by 4\/3 to, 
and that the approximation improves as to increases. 

We now consider the critical power of various vortex input profiles, and ask under what condition the critical power 
is close to the lower bound p cr (m). As in the vortex- free case, the only input profiles that can collapse with input 
power exactly equal to p cr (m) are ipo = XR m (Xr)e mL6 . We first calculate the critical power of the Laguerre-Gaussians 
profiles 

which are the vortex modes of the linear Schrodinger equation. To do that, we solve the NLS with the initial 
condition tpQ G and gradually increase c until, at Cth, the solution collapses. In this case, the critical power is close 
to p cr {m) for m = 1 but as to increases, the excess power above p cr (m) needed for collapse increases, see Tabic ILT1 
Similarly, for the sech input profile 

$» ech =cr 2 sech(r-5)e ime , 

the critical power is close to p cr (jri) only for m = 2, 3, 4, see Table ILT1 

To better understand these results, let us consider vortex profile of the form ipo = cf(r)e lm8 , where 

m=Q(p), P= r —^L : 

and Q{p) attains its maximum at p — 0. This ring profile is characterized by the ring width L and radius r max . As in 
the vortex- free case, the closer / is to a member of the one-parameter family XR m (Xr), the smaller the excess power 
above p cr (m) needed for collapse. By @, the family XR m (Xr) is characterized by 

radius/width = V3to. (3) 

Therefore, f(r) has to satisfy Q to "leading order" to be close to XR m (Xr). 

The Laguerre-Gaussian modes ipQ G are characterized by radius/width = -Jm/2. This ratio is close to © only 
for to ~ 1, explaining why the critical power of Laguerre-Gaussian modes is close to p cr (m) only for to = 1. Similarly, 
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TABLE II: Excess power above p C r(m) needed for collapse. 



Ill 


Input beams 




/ LG 




/ 7tl — sech 

Wo 


1 


0.65% 


20% 


0.13% 


2 


0.80% 


4.5% 


0.91% 


3 


7% 


1.9% 


0.71% 


4 


11% 


2.9% 


0.32% 


5 


14% 


9% 


0.17% 


6 


19% 


14% 


0.34% 



the sech profile ipQ ech is characterized by radius/width = 5. Since the radius/width of XR m (Xr) is equal to \/3m, this 
ratio is close to 5 for m = w 2.88, see equation ([3]). This explains why the threshold power of the sech profile tpQ ech 

is closest to p cr (m) for m = 3. As a final confirmation of this observation, we "fix" the sech profile ipff ch so that "it 
behaves like a XR m (Xr) profile", i.e., that it satisfies to leading order, as follows: 

^ra-sech = ^ (-^j SCch (r - VSrnj J™ 6 . (4) 

Indeed, the threshold power of the "modified" sech profile ^ is less than 1% above the critical power for m = 1, • • • , 6, 
see Table HU 

In Kruglov et al. estimated the critical power for vortex collapse to be equal to 

/(m) = 2 2m + 1 m!(m + l)! _ 
c (2m)! 

In they also estimated the critical power numerically for m = 1, 2, 3 and 4. These numerical results agree with our 
analytic calculation of p cr (m), but not with their own estimate ic , see Figure[T] To understand why this is the case, 

we note that the derivation of ii was based on the assumption that the collapsing vortex has a self-similar Laguerre- 
Gaussian profile. As noted before, the Laguerre-Gaussian modes are not a good approximation of the one-parameter 
family XR m (Xr), and as m increases this approximation becomes less and less accurate. In addition, the assumption 
that the solution undergoes an abcrrationless (adiabatic) self-similar collapse is known to lead to over-estimates of 
the critical power [5(. Indeed, even for Laguerre-Gaussian input beams, the critical power is closer to p cr (jn) than 
to I cm (m), see Figure [1] 

Most studies on optical vortices considered stationary vortices. Recently, there has been a growing interest in the 
dynamics of collapsing vortices. Bcrgc et al. showed that for vortices with input power P w ic , symmetry breaking 
noise causes the vortex ring to breaks into 2m + 1 filaments @. Vuong et al. generalized this result for vortices 
with power larger above ic" 7 ^ [Toj ] . We now show that these azimuthal instabilities can occur even for vortices with 

dimensionless power less than ic and even less than the lower bound p cr (m). To do that, we solve the NLS with 
the slightly elliptic Laguerre-Gaussian input profile 

= ^ G (Vz 2 + (1.05-y) 2 ) , (6) 

with m = 2 and with input power equal to |p C r("i = 2). Although the power of this vortex beam is below p cr (m), 
it breaks into two filaments which subsequently undergo collapse, see Figure [2J This effect of symmetry-breaking 
is very different from the case of peak-type non-vortex solutions, where deviations from radial symmetry increase 
the critical power for collapse This is because peak-type solutions collapse with the modulated Townes profile, 

(i.e., ip ~ T^TjR (xfj)) wnere L ~~ * a t the singularity) which is stable under azimuthal perturbations, as was 
demonstrated experimentally and numerically in [Tlj . and analytically in [l2T |. In contrast, vortices collapse with a 
ring profile, which breaks into a ring of filaments^nder azimuthal perturbations 0. Since these filaments do not 
collapse at the phase singularity point r = 0, each filament can collapse with the Townes profile, hence with the 
critical power p cr = p cr (m = 0) < p cr (m) (l3j . Note, that these filaments continue to rotate around r — 0, so that 
total hclicity is preserved. 
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FIG. 2: (color online) Solution of the NLS with input beam (j6j). Top: Levels set at z = 0,0.5 and z = 1 (from left to right) . 
Bottom: Surface plot at z = 1. 



Our results are also relevant for beams which are not linearly polarized. Let ip± be the amplitudes of the circular 
components e± = (x ± iy)/y/2. The equation for each circular component is 

i^- + A</>± + | [|^±| 2 + 2|^ T | 2 ] ^± = 0. 
In the case of a pure circular polarization (CP) state (?/>_ = 0), this equation reduces to 

i^+A^+ + H|V+lV+=0. (7) 

Since the Kerr effect is smaller by a factor 2/3 compared to the NLS (JTJ) for a linear polarization state, the critical 
power for collapse is larger by a factor of 3/2 [15j |. In particular, the lower bound for the critical power of a CP vortex 
beam -0+ = e lm Ao(r) is given by 

3 

Pcr P ( m ) = yPcr(m) PS 6^/3771. 

Similarly, consider the cases of radial polarization (RP) 

i> RP = A(r,t)[e ie e^+e- ie e+], 

and azimuthal polarization (AP) 

tP ap = iA(r, t)[e l0 e-~e- l9 e+\. 
Since = \ = \A\, the equation for each component is 

i^+AV±+2|V±|V±=0. 
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The Kerr effect is larger by a factor of 2, hence the critical power for collapse for each component is smaller by a 
factor of i, i.e., p cr (tp + ) = p cr (ip-) = hp cr (m = 1). In addition, the power of tp AP and ip RP is the sum of the power 
of ip+ and of ip- . Hence, 

P PP = Pcf = Pcr(i'+) +Pcr(lp~) = Pcr(m = 1) ~ 4.12p cr , 

in agreement with recent numerical simulations [l6|. 

In summary, we showed that the critical power for collapse of radially-symmetric vortex beams is typically a few 
percent above P cr (m) — 47r ^ „ 2 Pcr(™) where p cr (m) = R^rdr w 4v / 3ti. Deviations from radial-symmetry do 
not increase the critical power, but rather lead to disintegration into collapsing non-vortex filaments. 
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